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ABSTRACT 

We report the discovery of propeller-driven outflows in axisymmetric magnetohydrodynamic simu- 
lations of disk accretion to rapidly rotating magnetized stars. Matter outflows in a wide cone and is 
centrifugally ejected from the inner regions of the disk. Closer to the axis there is a strong, coUimated, 
magnetically dominated outflow of energy and angular momentum carried by the open magnetic field 
lines from the star. The "efficiency" of the propeller may be very high in the respect that most of the 
incoming disk matter is expelled from the system in winds. The star spins-down rapidly due to the 
magnetic interaction with the disk through closed field lines and with corona through open field lines. 
Diffusive and viscous interaction between magnetosphere and the disk are important: no outflows were 
observed for very small values of the diffusivity and viscosity. These simulation results are applicable to 
the early stages of evolution of classical T Tauri stars and to different stages of evolution of cataclysmic 
variables and neutron stars in binary systems. As an example, we have shown that young rapidly rotating 
magnetized CTTSs spin-down to their present slow rotation in less than 10^ years. 



1. INTRODUCTION 

Different accreting magnetized stars are expected to be 
in the propeller regime during their evolution. Exam- 
ples include accretion to fast rotating neutron stars (e.g., 
Davidson & Ostriker 1973; lUarionov & Sunyaev 1975; 
Stella, White, & Rosner 1986; Lipunov 1992; Treves, Colpi 
& Lipunov 1993; Cui 1997; Alpar 2001; Mori & Ruderman 
2003), white dwarfs in cataclysmic variables, and classical 
T Tauri stars (CTTSs) at the early stages of their evo- 
lution. The propeller regime is characterized by the fact 
that the azimuthal velocity of the star's outer magneto- 
sphere is larger than the Keplerian velocity of the disk at 
that distance. 

Different aspects of the propeller regime were inves- 
tigated analytically (Davies, Fabian & Pringle 1979; 
Lovelace, Romanova & Bisnovatyi-Kogan 1999; Ikhsanov 
2002; Rappaport, Fregeau, & Spruit 2004; Eksi, Hcrnquist, 
& Narayan 2005) and studied with computer simulations 
(Wang & Robertson 1985; Romanova et al. 2003; Ro- 
manova et al. 2004 -hereafter RUKL04). 

However, only relatively "weak" propellers were investi- 
gated earlier (RUKL04), in which a star spins-down, but 
no significant outfiows were observed. In this paper we re- 
port on axisymmetric (2.5D) simulations of "strong" pro- 
pellers, where significant part of the disk matter is re- 
directed to the propeller-driven outflows. We observed 
that the disk oscillates between "high" and "low" states 
and expels matter to conical outflows quasi-periodically. 
The quasi-periodic outbursts associated with the disk- 



magnetosphere interaction were predicted earlier by Aly 
& Kuijpers (1990) and observed in simulations by Good- 
son et al. (1997, 1999), Matt et al. (2002), Romanova 
et al. (2002 - hereafter - RUKL02), Kato et al. (2004), 
von Rekowski & Brandenburg (2004), RUKL04, Yelenina 
& Ustyugova (2005). However, none of the earlier simula- 
tions concentrated on the propeller stage, and only few os- 
cillation periods were obtained in earlier simulations. We 
report on modeling of propeller stage, where numerous os- 
cillations were observed. 



2. MODELING OF THE PROPELLER-DRIVEN OUTFLOWS 

We have done axisymmetric MHD simulations of the 
interaction of an accretion disk with magnetosphere of a 
rapidly rotating star. What is meant by rapid rotation is 
that the corotation radius of the star Vcr = (GM*/il^)^/^ 
is smaller than the magnetospheric radius which is de- 
termined by the balance between the pressure of the star's 
magnetic field and the ram pressure of the disk matter. 

The numerical model we use is similar to that of 
RUKL02 and RUKL04. Specifically, (1) a spherical coordi- 
nate system {t:,9,4>) is used to give high resolution near the 
dipole; (2) the complete set of MHD equations is solved to 
find the eight variables (p, Vr,ve,Vff,, Br, Be, B^, s) (with e 
the specific internal energy); (3) a Godunov-type numer- 
ical method is used; (4) special "quiescent" initial con- 
ditions were used so that we were able to observe slow 
viscous accretion from beginning of simulations (see de- 
tails in RUKL02). Compared to RUKL02, we now in- 
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Fig. 1. — Example of matter flow in the "high" state (at t/Po = 924). The color background shows the density distribu- 
tion which changes from red [p — 1.5) to blue {p = 0.0005). The solid black lines are magnetic field lines. The dark-blue 
lines/arrows are streamlines of the matter flow with the length of arrows proportional to velocity. The bold red line shows 
the surface where the matter pressure equals to magnetic pressure. 




Fig. 2. — The left-hand panels show the radial variation of the density p and azimuthal magnetic field i?0 above the equa- 
torial plane. The right-hand panels show the angular velocity 17 and the Keplerian angular velocity for comparison. 



elude the magnetic diffusivity. The viscosity and diffu- 
sivity are determined by turbulent fluctuations of the ve- 
locity and magnetic field (e.g., Bisnovatyi-Kogan & Ruz- 
maikin 1976) with both the kinematic viscosity vt and the 
magnetic diffusivity 7/t of the disk plasma described by 
a— coefficients as in the Shakura and Sunyaev model. That 
is, vt = ctvC^s/^K and rjt — adc^/riif , where VIk is the Ke- 



plerian angular velocity in the disk, Cs is the isothermal 
sound speed, and and ad are dimensionless coefficients 
< 1. In RUKL04 we investigated a range of small viscosi- 
ties and diffusivities, Ofv, ''^ 0.01—0.02 and found no sig- 
nificant matter outfiows. This paper investigates a wider 
range of a— parameters and finds substantial outflows for 
av ^ 0.1 and ^ 0.1 in the propeller regime. Enhanced 
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turbulence near the disk/magnetosphere boundary may 
arise because the radial gradient of the specific angular 
momentum is negative resulting in instability (Ustyugova 
et al. 2005). To model the diffusivity terms in the MHD 
equations, we used an implicit numerical scheme and the 
ICCG-method of solution of linear equations. 

2.1. Reference Units. The MHD equations were 
solved in dimensionless form so that the results can be 
applied to different systems. We take the reference mass 
Mq — Af*, a scale Rq = 2i?*, and a matter flux Mq which 
is close to the average matter flux through the disk. Then 
we derive a reference density po = ^J^o/voRq, a velocity 
vq — {GMo/Roy^'^ , a time-scale to = Ro/vq, and an an- 
gular velocity ^Iq = 1/to. The reference magnetic field 
is Bq = pqVq. Thus Bo = ^/pqVq. A reference magnetic 



moment is /io = ^o^o = vMqVqRq. The value of the 
magnetic moment used in our simulations /i* is typically 
10 times larger than /ig so that we introduce new reference 
variable /zqo = lO/iQ. The magnetic field at the surface of 
the star is i?, = Hoo/R^- The reference angular momen- 
tum fiux is A^o = MqVqRq. We measure time in units of 
the rotational period of a Keplerian disk at Rq, Pq = 27rto. 
We solve the MHD equations for the normalized variables, 
p — p/pa, V — v/vq, B = B/Bo, etc. and below show plots 
for normalized variables (with tilda's dropped). In para- 
graph 2.5 we show an example for CTTSs and millisecond 
pulsars in real units. 

2.2. Disk Oscillations and Outflov^rs. Here, we 
discuss results for a representative simulation run where 
outflows occurred. The parameters are /u* = hqq , ^* — 
(co-rotation radius rcr = Ro), c^v = 0.3, and ad = 0.2. 
The disk-magnetosphere interaction was found to be quasi- 
periodic. The system oscillates between a "high" state, 
where the inner radius of the disk is closest to the star, 
and a "low" state where the disk is at the largest distance 
from the star. Figure 1 shows an example of matter flow 
in the "high" state. The outflow is launched into a coni- 
cal shell of half-angle 6 ~ 45° — 60°. Many field lines are 
opened, and a major part of the matter flow to the wind 
is along the neutral line of the magnetic field. Analysis of 
the forces shows that the dominant force "pushing" mat- 
ter to the wind is the centrifugal force (Blandford & Payne 
1982). At smaller 9 there is a magnetically dominated out- 



flow of the much lower-density matter, which propagates 
with high-velocity along the open field lines of the star 
(see, e.g., Lovelace et al. 2002; Ustyugova et al. 2005). 

Our simulations show that: (1) matter accumulates in 
the disk and moves inward; (2) it comes close to the star 
and penetrates diffusively into the closed magnetosphere; 
(3) the disk matter acquires angular momentum from the 
rapidly rotating magnetosphere and is expelled as out- 
flows; (4) a small amount of matter accretes to the star 
through a funnel flow; (4) the disk is pushed outward by 
the rapidly rotating magnetosphere and the cycle repeats. 
Many field lines infiate and open during outfiow stage (see 
also Fendt & Elstner 2000) and some magnetic flux anni- 
hilates and may be a source of X-ray flares which are often 
observed in young stars (Feigelson & Montmerle 1999). 

Figure 2 shows the radial distribution of density p, an- 
gular velocity O, and azimuthal magnetic field B^j, in the 
"high" state (top panels) and "low" state (bottom pan- 
els). In the "high" state, the disk has a closest approach 
to the star of « 2.5, and the density in the inner disk 
is large. The magnetosphere rotates with super-Keplerian 
velocity. Azimuthal component of the field is small. In 
the "low" state, the disk is pushed outward, the density is 
lower, and the magnetic field lines are twisted up, forming 
modified expanded magnetosphere similar to the case of 
"weak" propellers (RUKL04). We conclude that no out- 
flows were observed in weak propellers (RUKL04) , because 
the diffusivity and spcciflc matter flux were too low. 

2.3. Efficiency of Propeller and Variability. We 
performed simulations for a range of stellar magnetic mo- 
ments p*, angular velocities $1*, and viscosity and diffu- 
sivity coefficients, and ad, and calculated the time- 
averaged matter fiuxes to the star (Ms) and to the wind, 
(M^), and the efficiency of the propeller. 
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The ratio TZ may be very large, 7?. ^ 50 — 100; that is, 
almost all of the matter coming inward in the disk may 
be ejected by the rapidly rotating magnetosphere. The 
dependence of 7?. on /i*, f2*, and is such that stronger 
ejection is observed in cases of stronger field, faster rota- 
tion and larger viscosity. However, the dependence on ad 
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Fig. 3. — The left-hand panel shows the matter fluxes to the wind and to the star for our reference case. The right-hand 
panel shows the quasi-periodic variations of the mass fluxes which we find for larger viscosity, in this case, = 0.6. 
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Fig. 4. — Color background and streamlines show fluxes of angular momentum carried by the matter (top panel) and by 
the field (bottom panel) in the "high" state (t/-Po = 924). Note that the scales are different, because Nf is very large 
near the surface of the star. However, the total integrated fluxes Nm and Nf have comparable values. 



has a turnover point at approximately 0.2. For ad <^ 0.2, 
TZ increases with ad, while for larger values it decreases 
(Ustyugova et al. 2005). This paper shows the dependen- 
cies only for ad ^ 0.2. 

It is important to note that no significant outfiows were 
observed when the diffusivity was relatively low, ad 0.1. 
This is because at low ad the penetration of the disk 
matter into the magnetosphere is not significant. Fur- 
ther, no outflows were observed at low viscosity. Analysis 
of the stresses show that the viscous stress is largest at 
the disk-magnetosphere boundary, and thus it adds to the 
"friction" and angular momentum transport from magne- 
tosphere to the disk. From other side, the matter fiux 
in the disk is proportional to the and at larger a^, 
disk penetrates to deeper, faster rotating layers of mag- 
netosphere. Both factors are important in generation of 
outflows, which appear at av ^ 0.1. Outflows were ob- 
served at a wide range of the magnetic Prandtl numbers, 
Pm = i^t/vt = ctv/am ~ 0.2 — 6 and thus do not require 
the dominance of viscosity or diffusivity. Instead, both a 
parameters should be larger than 0.1. Note, that the ob- 
served oscillations are completely determined by the pro- 
cesses at the disk-magnetosphere boundary. They are dif- 



ferent from the viscous instability oscillations of the disk 
(Kato 1978) which can not be investigated by our numer- 
ical model. 

The amplitude of the fluxes changes rapidly (see Fig- 
ure 3). There is a typical time-scale of variations Tqpo in 
each case which depends on main parameters fi*, ad, 
and av. For given ad and a^, the "quasi-period" increases 
with /i, and Q^. It varies in the range Tqpo = (5— 100)Po. 
We observed that for values of the diffusivity ad = 0-2, 
but relatively high viscosities ranging from 0.6 to 1, 
the oscillations become highly periodic. In one of of sam- 
ple runs a quasi-period changes from Tqpo = IOPq to 6.5Po 
(see Figure 3, right panel). Period changes due to the fact 
that the inner disk radius moved closer to the star. 

2.4. Angular Momentum Transport and 
Spinning-down. The angular momentum flux carried by 
the disk matter is re-directed by the rapidly rotating mag- 
netosphere to the outflows (top panel of Figure 4). Fur- 
thermore, there is a strong outflow of angular momentum 
(and energy) carried by the twisted open magnetic field 
lines from the star (Nf), the Poynting flux, and the closed 
field lines connecting the star and disk (bottom panel of 
Figure 4; see also Lovelace et al. 2002). The time aver- 
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aged total angular momentum flux from the star is (Ns) ~ 
-3.1(/i,//ioo)i-i(f^*/f^o)'-"K/0.2)"-«K/0.2)"-i. For 
our typical parameters the spin-down associated with 
open and closed field lines are comparable. However, at 
larger fi* and/or /i*, the outflow along the open flcld lines 
dominates (see related cases in Lovelace et al. 1995; Matt 
& Pudritz 2004), while at lower fi* or /u*, the situation 
reverses and a larger flux is associated with the closed field 
lines (like in Ghosh & Lamb 1979; RUKL02). 

The spin-down time-scale follows from equating the 
torque {Ns)No to I^,d^l^, /dt, where « lO^^gcm^ is the 
moment of inertia of the star. For the period of the star 
P* = 27r/^l*, we obtain: P*(f) = P*(0)(1 -|- t/tsd), where 
the spin-down time is 

2.5. Young CTTSs and Millisecond Pulsars. Our 

results can be directly applied to stars with relatively small 
magnetospheres, ^ (3— 10)i?*, for example, to CTTSs 
or to accreting millisecond pulsars. Thus, for a CTTS 
with a mass M* = CSM©, radius i?* = 2P0, an accre- 
tion rate Mq w 5 x 10~^ Mo/yr, and other typical pa- 
rameters Q!d = 0.2, Q!v = 0.2, and /i* = /ioo, wc obtain 
tsd ~ 5.8 X 10^ yr. We also derive dimensional values: 
/u* w 6.2 X 10^6 Gcm^, w 2.2 x 10^ G and P* w 1 day. 
Thus, if young CTTSs have strong magnetic flcld, then 
they spin-down rapidly to their presently observed slow 
rotation rate. 

Accreting millisecond pulsars have a different history of 
evolution: they spin-up from slow to fast rotation. How- 
ever, they may have episodes in the propeller regime. For a 
neutron star with mass = IAMq and similar accretion 
rate Mq ~ 5 x 10~^ M© /yr, we find tsd ~ 10^ yr. Taking 
R* = 10^ cm and /U* = /xqo, we find /x, « 7.0 x 10^'' Gcm^, 
P, « 7x 10^ G, P* « 1.3 ms. In this case tsd represents the 



time-scale of spin change for rapidly rotating millisecond 
pulsars. 

3. CONCLUSIONS 

In the propeller regime of disk accretion to a rapidly ro- 
tating star, we find from axisymmetric MHD simulations 
that the disk oscillates strongly and gives quasi-periodic 
outflows of matter to wide-angle (x w 45° — 60°) conical 
winds. At the same time there is strong field-dominated 
(or Foynting) outflow of energy, angular momentum and 
matter along the open field lines extending from the poles 
of the star. Matter outflows with high velocities and is 
magnetically driven. The outflows occur for conditions 
where the magnetic diffusivity and viscosity are signifi- 
cant, Qiv.d ^ 0.1. For smaller values of the diffusivity, the 
disk oscillates but no outflows are observed (RUKL04). 
The observed oscillations and outbursts are a robust re- 
sult, based on a numerous simulations at different param- 
eters with more than a 100 oscillation periods observed in 
many runs. The period of oscillations varies in different 
runs in the range Tqpo ^ (5 — 100)P*. It increases with 
/X* and ri*. We observed that the oscillations for rela- 
tively large become highly periodic with definite quasi- 
periods. More detailed analysis of these features will be 
reported later. A star spins-down rapidly due to both the 
disk-magnetosphere interaction and the angular momen- 
tum outfiow along the open field lines. The results are ap- 
plicable to young CTTSs, neutron stars, and cataclysmic 
variables. 
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